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Complete State Reconstruction of a Two-Mode Gaussian State via Local Operations 

and Classical Communication 

Gustavo Rigoliro and Marcos C. de 01iveir4H 

Instituto de Fisica "Gleb Wataghin" , Universidade Estadual de Campinas, 
Unicamp, 13083-970, Campinas, Sao Paulo, Brasil 

We propose a strictly local protocol completely equivalent to global quantum state reconstruction 
for a bipartite system. We show that the joint density matrix of an arbitrary two-mode Gaussian 
state, entangled or not, is obtained via local operations and classical communication only. In contrast 
to previous proposals, simultaneous homodyne measurements (HM) on both modes are replaced by 
local homodyne detections and a set of local projective measurements. 

PACS numbers: 03.67.-a, 03.67.Hk 



The feasibility of a quantum information task is related 
to the reduced or absence of non-local resources needed to 
its implementation and is an important asset for quan- 
tum communication purposes [l|, setting the limit for 
its widespread use. However, Quantum State Tomogra- 
phy (QST) |2j, y], a key tool in quantum information, 
is performed mostly through non-local operations. QST 
is a complete state reconstruction scheme implemented 
through a set of measurements over an ensemble of iden- 
tical quantum systems. For qubit systems it corresponds 
to the determination of all the Stokes parameters |4|] . For 
Gaussian continuous variable (CV) systems, as given by 
quantized electromagnetic field modes, it stands on a set 
of joint quadrature measurements, from which the joint 
density matrix p is reconstructed. Thus for Gaussian 
states, QST is equivalent to the measurement of global 
covariance matrices of the modes. For a two-mode Gaus- 
sian state most QST protocols to date either require si- 
multaneous HM on both modes [a, |y, |7| , with an exquisite 
control of both local oscillators (LO) phases, or require 
previous non-local operations on the modes to achieve a 
complete state reconstruction [8(. It is desirable, there- 
fore, the construction of a QST protocol that does not 
require any non-local operation and no phase-locking. In 
other words, a process which is operationally equivalent 
to QST, but without unnecessary non-local resources to 
its implementation. 

In this paper we show how one can reconstruct the 
whole density matrix p\ 2 of an arbitrary two-mode Gaus- 
sian state via local operations and classical communica- 
tion (LOCC) only. Since simultaneous HM of the two 
modes (a, l6| are not required, there is no need for con- 
strained control of the LO's phases, thus increasing the 
overall efficiency of the protocol, and also reducing the 
computational post-processing of data (See Ref. [9( for an 
interesting single homodyning alternative scheme). In- 
stead, a set of local parity and vacuum projections plus 
local squeezing are required. Our protocol is built basi- 
cally on three premises: (i) Alice and Bob can implement 
independent single mode local QST, certifying that they 
have a Gaussian state. Actually, after confirming (or be- 



ing informed previously) that one deals with a Gaussian 
state, only HM's of the variances of the modes will suf- 
fice, (ii) Both Alice and Bob are able to implement local 
squeezing and a local rotation on the quadratures of their 
modes. (Hi) Bob (or Alice) can make two types of local 
measurements: even/odd parity projections and vacuum 
projections of his (or her) mode. 

A bipartite two-mode Gaussian state pu is completely 
described [lfj, [ll| by its Gaussian characteristic func- 
tion C(z) = Tr(D(z)p 12 ) = e-^ z>Vz , where «t = 
(zi,Zi,z%, Z2) arc complex numbers, D(z) = e~ z Ev is 
the displacement operator, with E = disig(Z, Z), Z = 
diag(l,-l), and v = (vi, v 2 , v 3 , v 4 ) T = (01, a\, a 2 , a\) T 
the annihilation and creation operators of modes 1 and 2, 
respectively. T is the transposition, so that v is a column 
vector, and we have assumed all the first order moments 
to be null [121 ]. The covariance matrix V describing all 
the second order moments Vij = (— l) l+J (vii>J + v^Vi}/2 
is given by 
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Here Vi and V2 are the local covariance matrices of 
modes 1 and 2, respectively, giving the local properties of 
the two modes while C is the correlation between them. 
Finally, in addition to being positive semidefinite, V > 0, 
a physical Gaussian state must satisfy the generalized un- 
certainty principle, V + ^E > [id ]. 

The main goal of Alice and Bob is to obtain via LOCC 
the matrix V. Therefore, the first logical and trivial step 
consists in the measurement of Vi and V2 by Alice and 
Bob, respectively. These two covariance matrices are lo- 
cally obtained via any standard single mode HM tech- 
nique (or local QST). Up to now no classical commu- 
nication is needed and only after finishing this task Bob 
(Alice) informs Alice (Bob) of his (her) result. It is worth 
noting that we assume Alice and Bob have at their dis- 
posal a trustful source, in the sense that it produces as 
many as needed identical copies of the two-mode Gaus- 



sian state. 

The next non trivial step is the determination of 
C. To achieve such a goal, Alice and Bob need to 
work collaboratively [13[. First, on a subensemble of 
the copies, Bob implements parity measurements on 
his mode and informs Alice the respective outcomes 
for each copy, i.e., even parity (even number of pho- 
tons) or odd parity (odd number of photons). With 
this information Alice separates her copies in two dis- 
tinct groups, the even (e) and the odd (o) ones [l3(, 
as depicted in Fig. [TJ Alice's even group can be de- 
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FIG. 1: (Color online) Left: Alice separates her copies in two 
groups conditioned on an even (green) or odd (red) parity 
result obtained by Bob. Right: Alice selects copies corre- 
sponding to Bob's no-photon results. 

scribed by the non-normalized density matrix p\ = 
Tr 2 jP|p 12 P|} = £^=0 2 (2n|p 12 |2n) 2 , where P| =X X <g> 
S n =o |2^)2 2(2"|, T-\ the identity operator, and |n) 2 the 
n-th Fock state for mode 2. Using a similar notation, 
Alice's odd group is given as p\ = Tr 2 {P£ pi 2 P 2 } = 
S^Lo 2 (2" + l|pi2|2n + 1) 2 . But one can show that 
[3 cti = 2 v /dct(V 2 )( /5 f - pi) = Jdz lC < z ^c-^ r ^, 
where dzi = (l/7r)dRc(zi)dIm(zi) and Ti is the Schur 
complement |l_5| of V2 : 



Ti =V X - CV^C 1 = 



m mi 

Mi Vi 



(2) 



However, any one-mode Gaussian operator can be writ- 



ten as I dzie Zl *e 2 Zl lZl , being Ti its covariance ma- 
trix [lOJ I . Therefore, <j\ is a Gaussian operator whose co- 
variance matrix elements are r\\ = 2y / det(V 2 )((ajai} e — 
(a|ai) ) and p x = 2 v /dct(V 2 )((a 2 ) e - (a?) ), where 



(•) e = Tr(-pf) and (-) = Tr(-p°). Summing up, Ti can 
be obtained with the knowledge of det(V 2 ) and the sec- 
ond moments of p\ and p°, all of which determined via 



locc [lj 

S = (mi - 



Defining 7 = (?ii — 771) 



\m 2 \ 



and 



Ml) 



|m 2 | 2 ), Eq. © gives two indepen- 



dent equations, which alone cannot give m s and m c un- 
equivocally: 

(|m c | 2 + \m s \ 2 ) -2Re(m 2 m s m*), (3) 

7l 2 m c — 77l 2 77l s . (4) 



7 = n 2 (\m c 
S = 2n 2 m s m 



A unique solution though can be obtained if we con- 
sider an additional subensemble on which Bob performs 



another kind of projective measurement. The results of 
this measurement are communicated to Alice who build 
a local covariance matrix that is related to the origi- 
nal one through the Schur complement structure, sim- 
ilar to ^. In the present case we consider the simplest 
choice, i.e., Bob is able to perform a vacuum state pro- 
jection on his copies: photon-number measurements with 
no outcome. For each measurement, Bob informs Al- 
ice to which copies a no-photon result (p 2 —* |0} 2 2(0|) 
occurred. Alice, then, proceeds in a similar fashion 
as before but considering only the vacuum projected 
subensemble (right of Fig. [JJ), described by the density 
matrix p\ ac = Tr 2 (|0) 2 2 (0|pi 2 )/Tr 12 (|0} 2 2 (0|pi 2 ). One 
can show that [H| p\ ac = Jdzie z i Zai e~5 ZirilZl , where 



IIi = Vi C V 



b 



£1 V\ 



(5) 



with I the identity matrix of dimension two. Here, £1 = 
(a\ai) va c and v x = (a\) vac , where (-) vac = Tr(-p v 1 ac ). 
Explicitly, Eq. ([5]) gives us two more equations, 

a = (n 2 + -J (\m c \ 2 + \m s \ 2 ) -2Re(m 2 m s m* c ),(6) 

(3 = 2 f n 2 + - J m s m c - m 2 m 2 - m 2 m 2 , (7) 

in which a = (ni-£i)((?i 2 + l/2) 2 -|m 2 | 2 ) and f3 = (mi — 
vi)((n 2 + 1/2) 2 — |m2| 2 ). It is worth to stress that a and 
(3, as well as 7 and <5, are functions of parameters locally 
obtained by Alice and Bob. In order to solve Eqs. ([3]), ([!]), 
([rj|), and ([7]) for m c and m s we write rrij = \irij\e iej , where 
j = 1, 2, c, s and 9j is real. In this notation, our task is to 
determine 9 C , 9 S , \m c \, and \m s \ with the aid of T± and 
IIi, obtained by Alice via LOCC. The other quantities, 
rii, n 2 , |toi|, \m 2 \, 9\ 1 and 9 2 are easily obtained via local 
HM of modes 1 and 2. 

(i) Determination of 9 C and 9 S . Subtracting Eq. ((4]) 
from ([7]) we have f3 — S = m s m c , which gives 



\m s m c \ = \f3-5\, 

6 S + 6 C = Argifi-S), 



(8) 
(9) 



where Arg(z) is the phase of the complex number z. By 
the same token, subtracting Eq. ^ from © we have 



K| 2 + |m s | 2 = 2(a- 7 ). 



(10) 



Inserting Eqs. §E§ and (flU)) into © we get, a = (2n 2 + 
l)(a - 7) - 2|m 2 (/3 - 6)\ cos(9 2 +9 S - 9 C ). We could have 
used Eq. (J3]) as well. Solving, then, for 9 S — 9 C we obtain, 

9 S -9 C = cos^Ka^ - 7 (7i 2 + l/2))/|m 2 (/3 - 5)\] - 9 2 . 

(11) 
Eqs. j9|) and (fTTj) can be easily solved to give 9 C and 9 S , 
the phases of m c and m s . It is worth mentioning that 
Eq. (fTTj) is only valid when \(J3 — S)m 2 \ ^ 0. Later we 
show how to overcome this limitation. 



(ii) Determination of \m c \ and \m s \. From Eq. $W\ 



wc note that if wc had \m c 



\m s \ the problem would 



be solved. Manipulating the real and imaginary parts of 
Eq. we get 

K| 2 -KI 2 = \/3\ain(9 -e c -e s )/[\m 2 \sin(6 2 -e c +e s )}. 

(12) 
Here 9p is the phase of (3. Eqs. (fTO)) and (p~2|) can be di- 
rectly solved to give \m c \ and |t77 s |, the moduli of m c and 
m s . Eq. (fT2|) is only valid for 777-2 sin(02 — $c + 9 S )\ 7^ 0. 
Thus, all the covariance matrix elements can be locally 
reconstructed with a set of appropriate measurements 
and classical communication, establishing the following 
important connection to Gaussian QST: Global QST is 
completely equivalent to local covariance matrix HM, lo- 
cal parity and vacuum state projections, and classical 
communication. This is our central result and in the 
rest of this Letter we show how the necessary conditions 
|(/3 — 6)m 2 \ ^ and |m 2 sin(# 2 — C + 9 S )\ ^ can always 
be obtained by the addition of local squeezing [17J . 

(Hi) Overcoming |(/3 — 5)m,2\ = or \1ri2 sm(0 2 — C + 
6 S )\ = 0. To properly solve these problems we must 
know which quantity is zero. The simplest check is im- 
plemented when Bob reconstructs V2, which allows him 
to know if 7Ti2 = 0. Alice and Bob can also discover if 
777 c = 777 s = (implying /3 — 5 = m s m c = 0) by testing if 
Vi = Ti = III, since the absence of correlation (C = 0) 
between the modes cannot change what the parties mea- 
sure locally (See Eqs. © and ©). Also, if Vi ^ T ± or 
Vi ^ III Alice and Bob are sure that C /0 and the 
first non-trivial check sets in. They must discover if either 
?77 c = and m s =/= 0, or 777 c 7^ and tt7 s = 0, or both 777 c 7^ 
and 777 c 7^ 0. If cither m c or 777,5 is zero it is obvious that 



\h\ = |det(C)| = ||777 S |- 
one can show [13[ that 
and using Eq. (fTCT 



But 



M = v / dct ( v 2)dct(v 1 -r 1 ) 

we see that if |J3| = 2(a — 7) we 
know for sure that either 777 s or 777 c is zero. If we do 
not have an equality m c 7^ and m s 7^ 0. For our pur- 
poses, as we explain below, we do not need to know which 
quantity, m c or 777 s , is zero |18j |. Finally, to discover if 
sin(6> 2 — C + 6 S ) = we use Eq. pip and the phase of 
777 2 . Of course, Eq. (fTTj) is only valid if |t77. 2 (/3 — S)\ 7^ 0. 
Therefore, if \m 2 ((3 — S)\ = we first need to solve this 
problem in order to test if sin(6*2 — 6 C + S ) = 0. Since now 
we know which parameter is zero we are ready to show 
how Alice and Bob can overcome this situation allowing 
them to use Eqs. © to |12|) to obtain C. See Tab. [Qfor 
an overview of the strategies to solve these problems. 
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If 7712 = the most general solution [19( is achieved 
implementing a local symplectic transformation (local 



quadrature squeezing and rotation) on mode 2 [20( , S = 
diag(Ii, S2), where Ii is a 2 x 2 identity matrix acting on 
system 1 and S2 is given as 



tS2 cosh 7-2 sinhr2 
sinh r 2 e lS2 cosh r 2 



(13) 



being s 2 and r 2 real parameters. The new correlation 
matrix V is connected to V by V = SVS^ [20(, or 
equivalent^ for j = l,2,Vj = SjVjS^C = SiCSV 
Applying S to ([T]), the off-diagonal term of V 2 is 777,2 = 
e -2is 2m2 cosll 2 r2 + m * sinh 2 r2 + e -is 2n2 sm h(2r 2 ). Set- 
ting s 2 = and using that 7772 = we have 



777 2 = 772sinh(2r 2 ), 



(14) 



i.e., a new covariance matrix with 777,2 7^ 0. After this 
operation we can proceed with the original protocol to 
reconstruct V, which can be transformed back to give 
V = S- 1 VSt" 1 , withS- 1 =diag(li,S 2 - 1 ) and 



° 2 — 



e lS2 cosh r 2 
— sinh r 2 



— sinh 7*2 
i~ lS2 coshr 2 



(15) 



If either 7tt, s = or m c = 0, or equivalently /3 — 5 = 0, 
we can obtain a new matrix C = SiCS^2 where both 
parameters are not zero via a local squeezing operation 
alone. This leads to 



m a 

fhr 



2 m s cosh r 2 + m c sinh r 2 , (16) 

lS2 m c cosh ?'2 + 777 s sinh ?'2 . (17) 



Setting s 2 = in Eqs. (fTB")) and |T7]) we see that 777 s 
and 777 c arc combinations of m s and tt7 c . Therefore, if 
777 s = or 777 c = the new coefficients are necessary 
different from zero whenever we apply a local squeezing 
operation on mode 2. As anticipated, we do not need to 
know which quantity was originally zero. As before, after 
this local transformation we proceed with the original 
protocol obtaining V and then V. It is worth noting 
that when the two situations occur simultaneously, i.e. 
7772 = and 777,,. = or ?77 c = 0, the same local squeezing 
operation solves at once both problems, as can be seen 
in Eqs. CLl, dD, and HID. 

Lastly, after being sure that \m 2 m c m s \ ^ we can 
proceed to test if sin(#2 — C + S ) =0 using Eq. (fTTj) and 
the phase of 7772 , all quantities locally determined. In case 
of a positive result, there exist three possible solutions. 
The first one is valid when 777,2 7^ mi and is achieved 
reversing the roles of Alice and Bob in the protocol, as 
discussed above. The remaining two possibilities, and 
more general, is to locally and unitary transform mode 1 
or mode 2 before we implement the protocol, in the same 
fashion as before. Therefore, wc need to show that there 
exists at least one local unitary operation acting on mode 
1 or mode 2 that eliminates such a problem. 



Let us begin with mode 2. Applying the symplec- 
tic local transformation S2 we get, after assuming that 
sin(6» 2 -9 C + 6 S ) = 0, 

tan(6> 2 + ^-6»c) = 2^±sin(6»2-s 2 )sinh(2r 2 )/S ± , (18) 



t 



A± 



' s I - 
|2\ 



(19) 



|ra 2 |(|m c | 2 + \m s \ 2 ) =F 2n 2 \m c m 
±3\m 2 m c m s \ — n 2 (\m c \ 2 + |m s | 2 ) 
+ (±\m 2 m c m s \ + n 2 (\m c \ 2 + |to s | 2 )) cosh(4r 2 ) 
±2\m 2 m c m s \ cos(2# 2 — 2s 2 ) sinh (2r 2 ) 



-A T cos(6> 2 



s 2 )sinh(4r 2 ). 



(20) 



Here A± (B±) stand for the two possible values for 
the cosine, i.e., cos(# 2 — 8 C + 6 S ) = ±1, respectively. 
From Eqs. (fT5|) and (fH))) we see that a local squeez- 
ing alone (r 2 7^ and s 2 = 0) on mode 2 can make 
tan(0 2 + 6s - C ) 7^ if m 2 is not real (8 2 ^ 0). How- 
ever, whenever m 2 is real a rotation on the quadratures 
(s 2 7^ 0) is mandatory. There is one last loophole to fix, 
namely, the rare instances in which A+ = (note that 
A- is always different from zero). This is fixed by al- 
lowing the other party, in this case Alice, to implement 
a local squeezing on mode 1. As shown in Eqs. (|16[) and 
(fl~7|) this operation allows Alice to change at her will the 
phases of m s and m c without altering # 2 , solving com- 
pletely the last problem. By the way, this is the other 
possible solution for the sin(6> 2 — 8 C + 6 S ) =0 case, i.e., a 
local squeezing directly on mode 1. 

In summary, we showed a strictly local protocol in 
which a two-mode Gaussian state is completely recon- 
structed without relying on simultaneous HM or non- 
local resources. Actually, the only resources needed for 
this protocol are the ability to perform single mode HM, 
local parity and vacuum projective measurements, and 
classical communication. We also showed the complete 
equivalence of this local protocol to QST for Gaussian 
states. This equivalence is important for quantum com- 
munication purposes since now we can achieve the same 
goals of QST without non-local resources and simultane- 
ous HM. The set of local parity measurements required 
here, however, may restrict the implementation of the 
protocol, apart from instances where this measurement 
can be, in principle, performed [221 ]. Finally, this new 
protocol raises several interesting problems yet to be 
solved. Firstly, it is unknown if a similar local protocol 
can be devised for more than two modes and, secondly, if 
there exist other optimal sets of measurements, than par- 
ity and vacuum projections, allowing the complete state 
reconstruction of a two-mode (or many-mode) Gaussian 
(or non-Gaussian) state in a simpler local way. 

GR and MCO thank FAPESP and CNPq for funding. 
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